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Abstract. The recently introduced local moment approach (LMA) is extended to 
encompass single-particle dynamics and transport properties of the Anderson impurity 
model at finite-temperature, T. While applicable to arbitrary interaction strengths, 
primary emphasis is given to the strongly correlated Kondo regime (characterized by 
the T = Kondo scale wk)- In particular the resultant universal scaling behaviour 
of the single-particle spectrum D(w;T) = F(-^; ^-) within the LMA is obtained in 
closed form; leading to an analytical description of the thermal destruction of the 
Kondo resonance on all energy scales. Transport properties follow directly from a 
knowledge of D(uj; T). The T/wx-dependence of the resulting resistivity p(T), which is 
found to agree rather well with numerical renormalization group calculations, is shown 
to be asymptotically exact at high temperatures; to concur well with the Hamann 
approximation for the s-d model down to T/ojk ~ L and to cross over smoothly to 
the Fermi liquid form p{T) — p(0) oc — (T/wk) 2 in the low-temperature limit. The 
underlying approach, while naturally approximate, is moreover applicable to a broad 
range of quantum impurity and related models. 
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1. Introduction. 

The Anderson impurity model (AIM) [|l], reviewed fulsomely in 0, remains as topical 
as ever. While traditionally the archetype for understanding [0] magnetic impurities in 
metals, and heavy fermion systems when coherence effects are suppressed by disorder or 
alloying, study of it has recently acquired further impetus with the advent of quantum 
dots [0,0], which form directly tunable mesoscopic realizations of the model. 

In the strong coupling domain of large local Coulomb interaction (U), the 
intrinsic low-energy physics of the AIM is of course that of the Kondo effect 0. 
Manifest dynamically in the many-body Kondo resonance arising in the single-particle 
spectrum D(u), the effect is characterized by a single low-energy Kondo scale, uk- 
In consequence, D(u>) exhibits scaling in terms of u/uk alone, with no explicit 
dependence on bare material parameters. What then is the form of this universal 
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scaling spectrum? That the question itself is obvious does not vitiate the difficulties 
involved in answering it. Conventional perturbation theory cannot in general handle 
strong interactions, necessitating the development of new, non-perturbative theoretical 
approaches; the difficulties being particularly acute when dealing with dynamical and 
transport properties. And while Fermi liquid theory determines the low-frequency 
form D(u)—D(0) oc — (\u\/uj k ) 2 , such behaviour is limited to the lowest of frequencies 
|u;|/Vk <C 1, and otherwise provides no clues to the general form of the Kondo resonance. 

We have recently considered this problem || via the local moment approach 
(LMA) f|-|§. This non-perturbative method, naturally approximate but physically 
and technically quite simple, has been developed to handle a broad class of quantum 
impurity models, as well as lattice-based problems such as the Hubbard or periodic 
Anderson models within the local framework of dynamical mean-field theory |T0|, [TTJl . 
For the normal (metallic host AIM considered in ||, the LMA provides a simple 
analytical description of the Kondo scaling spectrum on all energy scales u>/u K . It shows 
in particular that while Fermi liquid behaviour is naturally recovered at sufficiently low 
frequencies, the scaling spectrum is entirely dominated for |u;|/u;k ^ 1 by long tails that 
exhibit a very slow logarithmic decay; in contrast to either the power-law (Doniach- 
Sunjic) behaviour previously thought to arise from empirical fits to quantum Monte 
Carlo [|12|,[T3| and numerical renormalization group (NRG) |I3L|l3J calculations, or the 
simple Lorentzian form || inferred e.g. from crude extrapolation of low- frequency Fermi 
liquid behaviour. The resultant theory, while naturally approximate, was found to 
give very good agreement for essentially all frequencies with recent NRG calculations 

E3- 

The above comments refer of course to zero temperature (T), and the obvious 
question is whether the LMA can be extended to finite-T. If so, then for the strong 
coupling Kondo regime in particular, the single-particle spectrum D = D(u;T) should 
scale universally in terms of u/uk and T/uk, with uk again the T = Kondo scale; 
and the approach then affords an analytical handle on the (tu-differential) thermal 
destruction of the Kondo resonance. In addition, and importantly, knowledge of 
the T-dependence of the single-particle spectrum enables direct access to transport 
properties @, notably the resistivity p(T). 

It is these issues we consider here, focussing as in || on the symmetric AIM, 
since it is this case whose low-energy physics reduces in strong coupling to the usual 
Kondo model (with exchange but not potential scattering). The relevant background is 
introduced in §2; and in §3 practical extension of the LMA to finite-T is considered. For 
finite interaction strengths U = U/tiAq (with Ao the one-electron hybridization), the 
LMA captures the spectra on all energy scales, including e.g. the non- universal Hubbard 
satellites. In §4 a brief discussion is thus first given of the thermal evolution of the 
spectrum on all energy scales; followed by consideration of spectral scaling as the Kondo 
limit is approached with increasing U, where a scaling form 7iAqD(uj; T) = F ^^; 
is correctly recovered. 

It is the scaling regime on which we naturally focus in the remainder of the paper, 
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developing an analytical description thereof in §5. The asymptotic behaviour of the 
scaling spectrum in particular, encompassing high-frequencies for all temperatures and 
vice versa, is obtained explicitly; and, as for the T = limit ||, is found to be largely 
independent of the details of the LMA. The temperature dependence of the resultant 
LMA resistivity, which is found to agree rather well with fmite-f/ NRG calculations ||16|| , 
is considered in §6; where we show it to be asymptotically exact for T/uk 3> 1 fT7|j . 
to agree well with the Hamann formula |L8[ for the resistivity of the Kondo/s-d model 
down to T/ujk ~ 1, and to cross over to characteristic low-temperature Fermi liquid 
form (T^j p(T) — p(0) oc — (T/o;k) 2 for T/u^ <C 1. Some closing remarks are given in 
§7. 



2. Background. 

We begin with requisite background material, starting with the familiar AIM 
Hamiltonian [0,0: 

k,a (7 ' k,a 

The first term describes the non-interacting host band with dispersion e*,; while the 
second refers to the correlated impurity with on-site interaction U, and site-energy 
ei = — T£ for the particle-hole (p-h) symmetric AIM (for which the impurity charge 
n i = So-(^io-) = 1 f° r & M U and T). The final term in equation fl2.ip is the host- 
impurity coupling. 

In considering single-particle dynamics at finite-T, we focus naturally on the 
retarded impurity Green function G{u;T) (<-> G(t;T) = — i6 ) (t)({cj cr (t), c] CT })), with 
the single-particle spectrum D(co;T) = — ^Im G(u;T). In the non-interacting limit 
U = 0, G(u;T) reduces trivially to g(u) = [u + — A(cj)] -1 , with u + = u> + i0 + . 
Here A(u;) = Ylik^ikl^ ~ e fe] _1 ( = A^ioS) ~~ iAi(co>)) is the one-electron hybridization 
function such that Ai(u) = it J2k ViL^fa ~~ e fc) = ^i(— u )\ an d the usual hybridization 
strength is defined by A = Ai(u = 0) oc phost(^ = 0) (with to = the Fermi 
level). From a knowledge of D(u; T) generally, and hence the transport time t(uj; T) oc 
[A -D(co>; impurity contributions to transport properties follow via the transport 

integrals 0,[L6],^(]]. The resistivity p(T) is given in particular by 



p(Ty- ] 

P(0). 



dc x - (JJx^Il ] (2 . 2) 



with f(uj;T) the Fermi function; it will be considered explicitly in §6. The thermal 
conductivity, Hall coefficient etc. are likewise expressible p], [16], ^C§ as appropriate 
moments of t 1 (u>] T) df(u; T)/duj. 

The essential strategy behind the LMA [|5|-[|] is threefold, (i) Local moments, 
regarded as the first effect of electron interactions, are introduced explicitly from the 
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outset |]J. Despite the stark deficiencies of static mean-field theory (MF) by itself 
(i.e. unrestricted Hartree-Fock), it may nevertheless be used as a starting point for 
a genuine many-body approach encompassing the correlated electron dynamics that 
are the essence of the Kondo effect, (ii) To this end the LMA invokes a two-self- 
energy description that is a natural consequence of the underlying two degenerate, 
broken symmetry MF states; introducing non-trivial dynamics into the associated self- 
energies £<j(co>; T) via coupling of single-particle excitations to low-energy transverse spin 
fluctuations, (iii) The final key element behind the LMA for T = 0, likewise discussed 
further below, is that of symmetry restoration: self-consistent restoration of the broken 
symmetry endemic at MF level. If symmetry can be restored (as is always the case for 
the metallic AIM considered here), then self-consistent imposition of it reveals the finite 
timescale r ~ h/uj m on which it occurs; the associated energy u m oc uk being the Kondo 
scale, whose origin within the LMA thus stems directly from symmetry restoration. 
The impurity Green function is expressed formally as ]5|-|| 

G(uj ] T) = \Y. G ^ T ) ( 2 - 3 «) 

where 

G a (uj- T) = [u+ - A(u) - t a (u + ; T)]- 1 (2.3b) 

(and o =| / I or +/—). The corresponding self-energies T> a (u> + ;T) are separated as 

t a (uj + ; T) = ~U\fi\ + T, a (uj + - T) (2.4) 

into (i) a purely static Fock contribution (with local moment that alone would 
survive at MF level (the trivial Hartree piece of = ^ is canceled precisely by 
€i = — y); and (ii) an tu-dependent contribution T: a (uj + ;T) containing in particular the 
spin dynamics that dominate the low-energy physics of the problem. The G a (uj; T) may 
be cast equivalent ly as 

G a (u- T) = [g a (uY l - S CT (cu + ; T)]- 1 (2.5) 

in terms of the MF propagator 

g a (u) = [u+ + ^u\n\-H"T\ (2-6) 

and with S CT = a functional of Q^, Q^. [Note that the Q a (oj) here refer to, and 

S CT and G a are thus constructed from, (either) one of the underlying degenerate MF 
states; and that the rotational invariance of the problem is correctly preserved by the 
'spin-sum' equation ( |2.3a| ) as shown directly in |J.] Separating the retarded self-energies 

as 

£ CT (^+; T) = E*(u; + ; T) - iS a (^+; T) (2.7) 

analyticity demands that E^cu 4 "; T) be positive semi-definite for all uj and T; while p-h 
symmetry implies 

E CT (^+; T) = E I _ a (-a;+; T) E*(u,+; T) = -S R CT (-^+; T), (2.8) 
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such that (— -Im G a (uj] T) =) D a (u; T) = D_ a (—uj] T) and hence D(u; T) = D(—u; T) 
for the impurity spectrum. 

For T = it is of course more natural/convenient to work with time-ordered self- 
energies, but these are simply related to their retarded counterparts. We denote the 
T = 0, t-ordered Green functions by G t (uj) and G^(lj) (for which equations (2.3) hold 
with uj + = uj + iO + sgn(co>)); and the corresponding self energies by £ CT (w), such that 
H a {u) = — i sgn(u;)E^(u;) with E*(a>) > V to. The t-ordered and T = 

retarded Green functions are then related simply by Re G a (uj; T = 0) = Re Gj(u) and 
Im G a {oj] T = 0) = sgn(a;)Im Gj(a;); and likewise: 



S^ + ;T = 0) = S^( 



Si(a; + ;r = 0) = S^ 



The LMA includes in £ CT (a>) the non-perturbative 
class of diagrams shown in figure 1, that embody 
dynamical coupling of single-particle excitations 
to low-energy transverse spin fluctuations. These 
capture the spin-flip scattering essential to 
describe the strong coupling Kondo regime for 
U = U/ttAq 3> 1; while in weak coupling 
they ensure the LMA is perturbatively exact 
to/including second order in U |J. Other classes 
of diagrams may also be included 0, but are 
of minor importance; retention of the dynamical 
spin-flip scattering processes is however essential, 
and it is these on which we focus. Figure 1 
translates to 
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Figure 1. Principal 
contribution to the LMA 
T, a (ui), see text. Wavy 
lines denote U. 



(2.10) 



where U~ aa (t) = i(T(S~ a (t)S°)) is the (t-ordered, T = 0) transverse spin polarization 
propagator (shown hatched in figure 1). It is given at the simplest level by an RPA-like 
particle-hole ladder sum in the transverse spin channel; viz 



n 



TT CTCT M 



[1 - U°U 



Ott — era i 



10 



(2.11) 



with °Il~ cra (uj) the bare p-h bubble, itself expressed in terms of the broken symmetry 
MF propagators {Q a } (solid lines in figure 1). 

The final step in the LMA for T = |§-|9j is self-consistent imposition of symmetry 



restoration (SR), embodied in the condition T,^(u> = 0) = Ej 



0) at the Fermi level 



uj = 0; and hence for the p-h symmetric AIM (using equations (2.4,8,9)) by 
(£f(0+;T = 0)^) E?(0) = i[/M. 



(2.12) 
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For given U, equation fl2.12| ) amounts in practice to a self-consistency equation for the 
local moment \fi\ ||. In physical terms, its consequences in general are threefold ||[7]]. 
(i) Satisfaction of the SR condition guarantees that the low-energy properties of the 
system amount to a renormalization of the non-interacting limit, i.e. that the system is 
a Fermi liquid with well defined quasiparticle behaviour, (ii) Most importantly, solution 
of equation ( j2.12| ) generates a low-energy spin-flip scale u m that sets the finite timescale 



t ~ h/u m for restoration of the broken symmetry inherent at crude MF level (stemming 
in effect from dynamical tunneling between the degenerate MF minima). Manifest in 
particular as a strong resonance in Imn + ~(a;) centred (by definition of c«J m ) on u = u m , 
this is the Kondo scale; with u m oc exp(— irU/8A ) in strong coupling, (iii) If the SR 
condition equation (|2.12j ) cannot be satisfied then a doubly degenerate local moment 



phase results @]: the spin-flip scale u m = (as characteristic of the locally degenerate 
state) reflecting the fact that the broken symmetry /degeneracy cannot be restored 
(r = oo). This situation does not arise in the metallic AIM considered here, where 
SR is satisfied for all finite U = U/ttAq. But it is the self-consistent possibility of 
such inherent in equation (|2.12|) that enables the LMA to access the quantum phase 
transition from a Fermi liquid to a local moment state in problems such as the soft-gap 
AIM [f7, 15[] where the non- Fermi liquid local moment state arises. 



Before turning to the LMA for finite-T we note (for use in §5 below) that the 
conventional single self-energy £(u; + ; T), defined by G(u; T) = [uj + — A{uj)— £(cj + ; T)] _1 , 
is readily obtained as a byproduct of the two- self-energy description intrinsic to the 
LMA. Using equations (2.3) it is given for arbitrary T by 



1 

2 


'E T (o; + ;T)+E ; ( W + ;r)* 


- g{u)E^ 


-;T)t i (u + ;T) 


1 - lg(u) 











S(w+;T) = — t p * : 1 (2.13) 

1 - \g(u) S T (c 

with g(uS) = [oj + — A(u;)] _1 the U = propagator. 
3. Finite-T LMA. 

To extend the LMA to finite temperature we first obtain, and then consider the rather 
transparent physical content of, the T = retarded self-energy S -(c ( j + ;T = 0). This 
follows, using equations (2.7,9), from its t-ordered counterpart H a (u); itself given 
explicitly by equation ( |2.10|) , where the transverse spin polarization propagator therein 
satisfies the Hilbert transform: 

J n u)\- z 

— oo 

Here, z = co> + i0 + sgn(a;) for the t-ordered polarization propagator T\r aa (oj) contained in 
equation ( |3.1| ) (while z = u+i0 + [= uj + ] for the corresponding T = retarded propagator 
U- aa (uj + ;T = 0)); and x~ aa (u; T = 0) = ImU-^iuj) (= sgn(cc>)Imrr ,T ' J (c<; + ; T = 0)) is the 
corresponding spectral density of transverse spin excitations, such that x~ acr { UJ ] T = 0) > 
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and x~ aa (u;T = 0) = x a ~ a (-<*>; T = 0) [§. Using equation in equation (|2J0| ) 
leads straightforwardly to 

oo 

= u 2 J ^[eMg+^cu-^) + e{-u 1 )gz <T (u-cj 1 )] x^^-t^) (3.2) 



where 

y c<j — u>i ± iU + 

—00 

denote the one-sided transforms of the (t-ordered) MF propagator in terms of its spectral 
density D®(u); and 8(x) is the unit step function. The T = retarded self-energy then 
follows directly using equation (|2.9|), viz 



00 00 
Sa („+;T = 0) = f/ 2 |^ J du 2 



[0(wi)0(w 2 ) +9{-lo x )8{-u 2 )\ ' v— (^ i; T = 0) (3.4) 

a; — cJi — uj 2 + iU + 

as sought. 

The natural extension of the LMA to finite-T may now be deduced by considering 
the physical processes to which the two terms in equation ( ^.4[ ) for T = correspond. 
We focus explicitly on £| in the following; with f(u;T) = [exp(Pu) + the Fermi 
function, such that f(u; 0) = 9(— to). Consider the first term in equation ( ^.4[ ) for 
S^(ci; + ; T = 0), in 6{uji)6{uj2) (noting that x + ~( UJ i'iT) for uj\ > corresponds to first 
flipping the impurity spin from f to |). Physically, this corresponds to processes in 
which (for Ej_) a |-spin electron is first added to an |-spin occupied impurity, and 
the originally present |-spin then hops off the impurity/site into the conduction band. 
The latter generates an on-site |— >[ spin-flip, which is a hard core boson in the sense 
that at most one spin-flip can be created on the impurity; and the (T = 0) probability 
with which the |-spin can be added to/hop into (empty) conduction band states is 
9(u)2) = 1 — f(^2',T = 0). Appropriate extension of this contribution to finite-T is then 
achieved simply by replacing 6{uii)6{ijj2) with 6'(u;i)[l — f{uji\T)\ in equation ( |3.4| ) for 
Sj_; with 1 — f(tU2',T) now the finite-T probability with which the original f spin can 
hop into conduction band states, and to which process creation of the on-site spin-flip 
is again slaved. 

A directly analogous situation pertains to the second term in equation (3.4) for 



£j_, in 8(— uji)9{— LV2) (noting that x +_ (c ( j 1 ;T) for uj\ < corresponds to first flipping 
the impurity spin from | to f). This corresponds physically to processes in which (for 
an originally present J,-spin on the impurity is first removed, and an f-spin electron 
then hops from the conduction band onto the impurity. The latter process generates 
the on-site {— >f spin-flip; and the (T = 0) probability with which the |-spin electron 
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may be removed from the (occupied) conduction band states is 9(—u 2 ) = f(u 2 ;T = 0). 
Extension of this contribution to finite-T is then achieved by replacing $( — u!\)6( — uj 2 ) 
with 9(—uJi)f(u) 2 ;T) in equation ( |3.4|) for Ej_; with f(u 2 ;T) the finite-T probability 
with which an j-spin electron can be removed from the conduction band to hop onto 
the impurity, thereby generating the on-site spin-flip. 

The finite-T LMA we consider is thus given from the above by 



DC 



£^+;T) = f/ 2 / ^ / d^ X + -K;T) D *M) ( ) (3 5a) 

J IT J UJ — U)\ — U) 2 + lU" 1 " 

—00 —00 

with 

g(un u 2 ) = - f(u 2 ; T)] + 0(-Wi)/(u; 2 ; T). (3.5b) 

In equation ( |3.5fl| ) we have also naturally included the T-dependence of T) > 0; 

where again x +- (u;;T) = sgn(co>)Im IT (cj + ;T), with II + ~(cg> + ;T) now the finite- 
T polarization propagator (for which the Hilbert transform equation ( |3.1| ) is again 
applicable, with z = uj + ). Ef(co> + ;T) likewise follows from equation (3.5) by inverting 
all spins (T^l? + ^ — ) or, in the symmetric case of interest, from p-h symmetry 
equation (|2.8|) . 

We note that equation (3.5) is not what arises if one attempts to obtain 
the appropriate LMA £ CT (u; + ;T) via a conventional analytical continuation of the 
corresponding imaginary time H a (iu n ). In that case it is readily shown that the 
resultant E|(c<j + ;T) is indeed of form equation (|3.5c| ), but with g(ui;u 2 ) replaced by 



g(oJi;oJ 2 ) = -sgn(a;i)[/ B (-^i;T) + f(u 2 ;T)} where / B (o;) = [exp(f3u) - 1] 1 is the 
Bose function. This corresponds physically to treating the impurity spin-flips as free 
bosons, with thermal statistics that are entirely divorced from those of the fermions. In 
generating multiple thermally created spin-flips, it thus in effect violates the hard core 
condition for the bosonic spin-flips - whose probability of creation is dictated by the (T- 
dependent) probability with which fermions can hop to/from the impurity from/to the 
conduction band, as embodied in equation (3.5). The hard core constraint can however 
be recovered simply by replacing the free Bose function /b(— ^i',T) in g(ui,u 2 ) by its 
T = limit /b(— wi;0) = —6{uj\): in that case g(uji]uj 2 ) reduces precisely to g(ui,u 2 ), 
and equation (3.5) for Ej(cj + ; T) is recovered. We do not know what additional diagrams 
may be required to recover hard core behaviour within the framework of a conventional 
analytical continuation; but believe that the physical arguments given above, together 
with the results of the following sections, attest to the essential validity of equation (3.5) 
as used in practice for the LMA H a (uj + ; T). 

The finite-T LMA is readily implemented. The retarded G(uj; T) (and hence 
D(uj;T) = — ^Im G(u;T)) is given from equation (2.3), with the self-energies 
Eo-(o; + ;T) from equation (|2.4|) and Y la (uj + ;T) from equation (3.5). The finite-T 
retarded polarization propagator entering equation ( |3.5o| ) (such that x a ~ a {oj\ T) = 



sgn(u;)Imn o " a (w + ; T)) is calculated in practice for finite U at the level of the RPA-like 
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p-h ladder sum, obtained by straightforward analytical continuation of the imaginary 
time II + -(iw m ) = °n + -(icj m )[l - U°U + - (i^)]- 1 (cf equation pH|)). Finally, it is also 



straightforward to include a T-dependence for the local moment = |/x(T) | entering the 
MF propagators (equation (|2.6|)) and the LMA self-energies (equation (|2.4j)). This may 
be encompassed via \n{T) \ = |/i(0)| +5\fi{T)\, with |//(0)| the T = moment required to 
satisfy the symmetry restoration condition equation ( |2.12| ); and with 5\fi(T) \ calculated 
in practice at MF level, viz %(T)| = duo [D°(uo) - D®(oj)][f(u;T) - 9(-u)] with 
D° a (uj) = Aqit^/Kuj + ff/| / u(0)|) 2 + Ajj] the MF spectral density. We find however 
that the resultant T-dependence is negligible for essentially all U (provided one is not 
concerned with temperatures on the order of U), and we thus omit it from the results 
shown explicitly in §'s 4ff. 

4. Results. 

The primary interest in single-particle dynamics resides of course in thermal destruction 
of the low-energy Kondo resonance, particularly in the strong coupling Kondo/spin- 
fluctuation regime where the AIM maps onto the Kondo or s-d model. That is our main 
focus, and is pursued in the following sections. First however, we consider the issue of 
spectral scaling as the Kondo limit is approached with increasing U; preceded by a brief 
discussion of the thermal evolution of the spectrum on all energy scales. The following 
results, obtained from the finite-T LMA outlined above, refer explicitly to the usual 
wide-band AIM for which Ai(u) = A V to. Here and throughout, the Kondo scale 
c<jk is defined as the half-width at half maximum of the T = single-particle spectrum 
D(u;0); and uk oc u> m with the spin-flip scale u m determined (as above) via 

symmetry restoration, and exponentially small in strong coupling. 

Figure 2 shows the resultant spectrum, tiA q D(uj] T) vs u/A ; for a fixed interaction 
strength U = U/hAq = 4, and for temperatures T/u K = 0.1, 10, 40 and 160, the intent 
being an all scales overview of spectral evolution. The dominant effect of increasing 
temperature is naturally to erode the many-body Kondo resonance (as later considered 
in detail). This process occurs initially 'on the spot', with no effect upon spectral 



features on the non- universal energy scales to ~ A or U [|,16]; as evident in figure 2 
for T/ujz = 0.1 and 10, where e.g. the high-energy Hubbard satellites retain their T = 
form (being centred to high accuracy for the wide-band AIM on = ^ (l + ^f) U 
). By contrast, for non-universal temperature scales on the order T ~ A , spectral 
redistribution occurs on all energy scales including the high-energy satellites. This 
effect, likewise known to arise from NRG calculations [0, is evident in figure 2 for 
T/uj k = 40 and 160 (T/A ~ \ and 1 respectively). 

We turn now to the question of spectral scaling. In the strong coupling spin- 
fluctuation regime, the low-energy physics of the AIM depends solely upon the Kondo 
scale; and the T = Kondo resonance hence exhibits universal scaling in terms of uj/ujk 
alone, with no explicit dependence on the bare material parameters U or A . The 
LMA for T = leads correctly to such scaling behaviour with increasing U || |§ , and 
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Figure 2. vrA D(w;T) vs w/A for U = 4 and T/u K = 0.1 (solid line), 10 (long 
dash), 40 (short dash) and 160 (dotted line). 

the resultant scaling spectrum gives very good agreement with NRG calculations ]5|,|H)|]. 
Scaling behaviour should likewise arise in strong coupling for 0, but with the scaling 
spectrum nA D(ui; T) = F{oj /uk', T/uk) now dependent upon both u/uk and T/ujk- 

That such scaling indeed arises within the LMA upon progressively increasing U is 
shown in figure 3 where, for fixed T/ujk = 1, ttAqD(uj;T) is shown for three different 
interaction strengths, U — 4, 6, 8. The inset shows the spectrum on an absolute 
scale vs uj/Aq, illustrating the exponential narrowing of the Kondo resonance with 
increasing U (uk oc exp(— 7r[//8A )). The main figure by contrast shows 7rAo-D(co>; T) 
vs u/ujk- The spectra are indeed seen to approach asymptotically that of the Kondo 
scaling limit (obtained analytically in §5, and also shown in figure 3); which behaviour 
is reached in practice for the U's considered in figure 3. The scaling illustrated above 
arises generically, and figure 4 shows the resultant LMA scaling spectra irA Q D(uj] T) vs 
uj/uk for a range of different temperatures T/uk- 

Before proceeding to a detailed analysis of the scaling spectra, we comment briefly 
on the T-dependence of the transverse spin polarization propagator fI H (uj + ; T), whose 
spectral density enters the LMA self-energy S -(c<j + ;T) (equation ( |3.5q| )). Its 

T-dependence, arising from that of the °U + ~(uj + ; T) bubble, is readily shown to occur 
on the scale T ~ Ao and is non-negligible only for non-universal temperatures of this 
order (where it contributes e.g. to the 'spectral redistribution' discussed above in regard 
to figure 2). The Kondo regime by contrast corresponds formally to finite T/u>k in the 
limit uj k oc exp(— tcU /8A ) — > 0. Here the T-dependence of x + ~(y;T) is thus entirely 
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Figure 3. ttA D(lj;T) vs w/wk for fixed T/ojk = 1, and [7 = 4 (dotted line), 6 (short 
dash) and 8 (long dash). The approach to the Kondo limit scaling spectrum (solid 
line, see §5) is evident. Inset: Corresponding spectra on an absolute scale, vs w/A . 
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Figure 4. LMA(RPA) scaling spectra: 7rA L>(w;T) vs u)/w k for T/w K = 0.1 (solid 
line), 0.5 (long dash), 1 (short dash), 5 (point-dash) and 10 (dotted). LMA(RPA) 
denotes that x H ( w ) is given explicitly by the RPA-likc ladder sum (see §5.2). 
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negligible (so from now on is referred to solely as x + i 00 ))] an d the T-dependence of 
the LMA scaling spectrum is then controlled exclusively by g(u\\L02) (equation ( |3.56| )) 
entering the LMA S a (cu + ;T). 



5. Scaling Spectrum. 

Our aim now is to obtain analytically the u- and T-dependence of the strong coupling 
Kondo scaling spectrum arising from the LMA; and, as in || where the T = limit was 
considered, we do so with only minimal assumptions about the form of the transverse 
spin spectrum x + ~{ UJ )- In the following analysis the Kondo scale will appear in several 
guises, viz u m , u' m , wk and A Z (with Z the T = quasiparticle weight). These are all 
however proportional to one another, and thus equivalent; reflecting the fact that in the 
Kondo limit the problem is characterized by a single low-energy scale. 

Within the LMA the 'primary' manifestation of the Kondo scale is u m , the spin- 
flip scale determined via symmetry restoration (see §2 and below), and embodied in a 
strong resonance in x + ~{ UJ ) centred on u = ou m . And the scaling spectrum is obtained 
by considering finite u = oj/uj m and T = T/u m in the limit u> m oc exp(— ttU/8A ) — ► 0, 
which projects out the non-universal temperature and energy scales that are irrelevant 
to the scaling regime. Hence, referring to equation (2.3), the 'bare' u = uuj m = may 
be neglected; and likewise A(a>) = A(d;u; m ) reduces to A(0) = — iA (which applies to 
any metallic host, so that the following analysis is not confined to the wide-band AIM). 
The scaling spectrum is thus given from equation (2.3) by 

2 ^ (AJ'EJK; T)f + (1 + Ao-^K; T)f 

where Y> a (uj + ]T) = —o~^ + S cr (u; + ;T) in the Kondo limit where the local moment 
saturates (|/i| — > 1). The spectrum is naturally determined solely by the scaling 
behaviour of the self-energies, and the LMA £j(u; + ;T) is given explicitly by (cf 
equation (3.5)): 



oo oo 

^{uj + ] T) = U 2 J J du 2 



X + 'M - D ° iM ,. n+ [^(o; 1 )/(a; 2 ;r) + 0(-a; 1 )(l-/(a; 2 ;r))](5.2) 

In the strong coupling limit the transverse spin spectrum x H (u) has the following 
functional form 

-X + -{u) = —F{Co)d{Co), (5.3a) 
naturally scaling in terms of Co = u/uj m ; and 

oo oo 

J ^ x + -(^) = l = Ajdy F(y) (5.3b) 
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which reflects the saturation of the local moment and complete suppression of double 
occupancy in the Kondo limit. The function F(u) is distributed around and centred 
upon qj = 1 (by definition of u> m ), and F(u) ~ Cj as Cj — > 0. The above behaviour 
is readily shown to arise explicitly |6j with x + ^{ u ) given via the p-h ladder sum 
equation ( |2.11| ), for which F(u) has the functional form 

F(Q) = -. (5.4) 

K J l-2aCj + Cj 2 K J 

The specific form of F(Cj) is not however required in the following analysis (we shall 
need it only in §5.2); and the important asymptotics of the scaling spectrum considered 
in §5.1 are in fact independent of it. 

From equations (5.2,3), noting trivially that f(u;T) = f(Cj;T), it follows directly 

that 

oo 

V]{u+-T) = itU 2 A J dy F(y)f(y + uj;f)Dl(uj m [y + uj)) (5.5) 
o 

where D9(uj m [y + uj]) = -D|(0) since we consider finite Cj with uj m — > 0. But in the strong 
coupling limit, -D|(0) = — ^lm(^(0) is given (see equation ( |2.6| )) by 



irU 2 D° L (0) = 4A (5.6) 



and hence: 



oc 

Ao'^T^I T) = AAj dy F(y)f(y + Co; f). (5.7) 
o 

As required in the Kondo limit this scales solely in terms of uj and T, with no C/- 
dependence (A being a pure number 0(1) determined via equation ( |5.3£ )). 
From equations (5.2,3), E?(u; + ;T) is likewise given by 

oo ( oo "\ 

Ef (uj + ; T) = AU 2 [ dy F(y) \ [ du 2 D ^ f(uj 2 ; T) \ (5.8) 

J J a; m [cu + y\-uj 2 

^-oo J 

(with a principle value implicit from now on). For T = 0, the c^-integration leads to 
Re G7(uj m [uj + y]), with Re Q7 {uj) logarithmically divergent as u — » 0. To handle this 
at finite-T we proceed as in p[. The lower limit of the ^-integration in equation 



is replaced by a UV-cutoff, —A, where A = min[^,D] (with D the host bandwidth); 
its precise value being immaterial in the following. Rescaling uj 2 to uj 2 = uj 2 juj m in 
equation (|5.8j) then leads to 



VSf(« + ;r) = [dy F(y) [dCj 2 J^f} (0.9) 

7f J J uj 2 -[uj + y\ 

o —A 
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where equation ( |5.6| ) has again been used; and A = \/uj m is thus defined, such that 
A —>■ oo in the Kondo limit. 

For T = 0, equation ( [5~9|) reduces to 



A^Eftw+j 0) = - In 

7T 



" A " 




_U) m _ 










(5.10) 



(where equation fl5.3fcp is used). From this the explicit [/-dependence of the strong 
coupling Kondo scale u m follows directly from symmetry restoration (equation ( f2.12j )), 
viz S?(0; 0) = \XJ\ namely 



where 

uj' ni = Aexp(-7r[//8A ) 
and c a constant 0(1) given by 



c = exp 



oo 

-A J dy F{y) my 



(5.11a) 



(5.116) 



(5.11c) 



The prefactor to the exponent of u m or u' m is of course approximate (reflecting the 
UV-cutoff used above), but the exponent itself is exact. And from equation (|5.10|) , 
Ef {lu+; 0) = Ef(u+; 0) - Ef (0; 0) at T = is given by: 



AA 



A o - 1 Ef(cu + ;0) = / dyF(y)\n 

7T 



1 + ^ 



(5.12) 



For the scaling spectrum D(u>;T) (equation (|5.1|)) we finally require E^(u + ;T) = 
-| + £f(u;+;T), i.e. 

Ef(oo + ; T) = [Ef T) - Ef (uj + ; 0)] + Ef 0) (5.13a) 

with E^(u) + ;0) given above. And using equation (|5.9| ), a straightforward calculation 
yields 



4A 



A, 1 [Ef(u + ;T) - Ef(u + ;0)} = -— / dy F(y)H 



\y + w\ 
f 



(5.136) 



where 



H{x) 



duj P 



u) + x u) — x J exp(u;) + 1 



(5.14) 



and the limit A — > oo has been taken with impunity. As for A 1 Ef (tu + ; T) 
(equation (|5.7| )), equations (5.12,13) show that Aq 1 E?(o; + ; T) also scales solely in terms 
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of Cj and T with no U -dependence, as required from equation ( |5.1|) for the spectrum to 
scale thus. 

Equations (|5~7j) and (5.12,13), together with p-h symmetry (equation ( f2.8|) ), provide 
the basic results from which the Cj- and T-dependence of the scaling spectrum may be 
determined; we analyze them in the following sections. 

Before proceeding we note that the function H(x) (equation ( |5.14j )) is given in 
closed form by 



H(x) = - In 



\x\ 
2^ 



(5.15) 



with ip(z) = dlnT(z)/dz the digamma function; and is readily shown to be of form 



H(x) 



In 



\X\ TT 



1 2\ 7 



h{x) 



(5.16a) 



where h(x) is non-singular. A useful approximation to h(x) is 



h(x) 



7 



1 + §7 [fl 



■ 1 2 



(5.166) 



where 7 = C + In 2 and C (= 0.577...) is Euler's constant. As required in the following 
section, this ensures that the leading asymptotics of H(x) as x — > and x — > 00 are 
captured exactly; being given by 



and 



H(x 



H(x) 



-In 



6 



\x\ 

7T 

'7T 

X 



- 7 + 0{x 2 



0(x~ 



(5.17a) 



(5.176) 



(which is essentially a Sommerfeld expansion). 



5.1. Spectral asymptotics. 

We now consider the asymptotic behaviour of the LMA scaling spectrum and, relatedly, 
the corresponding single self-energy £(u; + ; T). This encompasses both the spectral tails, 
\uj\ ^> 1 for any T, as well as high temperatures, T ^> |1 + Q\ for any uj\ and is not 
dependent on details of the function F{uj) (equation (5.3)) that determines the spectral 
density of transverse spin excitations. In practice the specific form of F(u) is required 
only for low u and T, as considered in §5.2. 

From equation ( |5.7|) , using equation (|5.36|) and recalling that F{Cj) is distributed 
around and centered on Co — 1, Aq 1 Y^{uj + ] T) is given for \Cj\ 1 and any T by 

At 1 V\(u + ] T)~4:f{u> ] f). (5.18) 
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This holds also for T ^> | and all uj. Turning now to equation (5.13) for T), 
note first from equation (|5.12|) that for T = 0, Aq 1 S| l (ci; + ; 0) ~ — ^ J °° dy F(y) In \u>/y\ 
for |u)| ^> 1; and hence, using equations (5.11a,c) and (|5.3<j) , 

4 



with 



V2f(o> + ;0) 



= uj/uj' m 



~ In I a;' I 

7T 



(5.19) 



(5.20a) 



thus defined. Equation ( 5 . 1 3 6|) for \uj\ ^> 1 likewise reduces (again using equation ( |5.36| )) 
to A^pf^T) - Xf(u;+;0)] ~ -J#(|£|/T) = -ftf(|u/|/n where 



T' = T/J m 



(5.206) 



is correspondingly defined. Combined with equation (5.19), and using equation ( 5.16o| ) 
for H(x), equation (5.13) for T^(u + ; T) thus reduces for \uj\ ^> 1 to 



A - x Sf (w+; T) 



~ ( In 

7T 



a; 



/|2\ * 



) 



h 



(5.21) 



which is readily shown to hold also for T ^> |1 + u\ and all Co. 

The asymptotic behaviour of the Kondo limit scaling spectrum then follows from 
equation (|5.1|), using equation (|5.18|) (and p-h symmetry), as 



ttA D(uj;T) ~ (5.22) 

if [l + 4/(^;TQ] [l + 4/(-^;TQ] 

2 1 [A - 1 Sf(^;r)]2 + [l + 4/(^T')] 2 [Aq 1 E?(o;+; T)] 2 + [1+4/(-u/; T')} 2 



with A 1 E 1 | l (ti; + ; T) from equation (|5.21j) . Three points should be noted here: 
(i) For \uj'\ > max(l,T') (where f(u';T) ~ 0(-c</))» equation (ggg) for Zf(tu + ;T) 



reduces (using equation (|5.16fl )) to its T = limit, equation (|5.19|) ; hence the spectrum 
for \u'\ ^> max(l,T') reduces to 

irA D(co; T) ~ - \ ?— = + ^ I . (5.23) 

2\[iln M ] 2 + l [llnM] 2 + 25j 

This is also the result of || for the T = spectral 'tails', there shown to dominate 
the T = scaling spectrum (down to uo ~ c^k, the crossover to Fermi liquid behaviour 
occurring only on the lowest energy scales \u)'\ <C 1); and to be quantitatively accurate in 
comparison to NRG calculations for \u'\ > 5. The arguments leading to equation ( |5.23j ) 
show that the same behaviour arises also at finite-T, for sufficiently high frequencies 
|u/| 3> max(l,T'). This 'common tail' behaviour is indeed seen in the scaling spectra 
shown in figure 4; and is evident also in recent finite-T NRG calculations for the Kondo 
model itself ( figure 1 therein). 
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(ii) For X" ^> |1 + uj'\ by contrast (where f(u>'] T') ~ £), equation ( |5.21| ) for Y^{uj + ; X) 
reduces (again using equation ( |5.166| )) to 

A x tf (uj + ; X) ~ -- ln[eT'] (5.24) 

7T 

where e = |exp(— C); and the spectrum is thus given by 

nA D(cu;T) T J- T - r -. (5.25) 

V ; 3 + ^ln 2 [eX'] V ' 

For X' ^> |1 + uj'\ the cu-dependence of D(u; X) is thus essentially irrelevant; as is seen 
in figure 4 for e.g. T/uk = 10, where the spectrum flattens out at low frequencies. 
Equation ( |5.25 ) gives in particular the asymptotic high-temperature behaviour of the 
Fermi level spectrum XJ(0;X), which will be considered further in §5.2 

(iii) For X' 3> 1, the leading asymptotic behaviour of D(u; X) is clearly dominated by the 
logarithmic growth of X) (since /(u/; X") < 1); and is given from equation (|5.22j ) 
by 



3tt 2 



16 



In 2 



V^vrX') 2 + |u/| 2 e" 



(5.26) 



which form will prove important in §6 where the resistivity p(T) is considered. 

Finally, the corresponding asymptotics of the conventional single self-energy 



£(cu+;X) 



X) — iE 1 (uj + ;T) may also be obtained. S(cj + ;X) is related 



generally to the LMA self-energies T: a (uo + ;T) by equation (|2.13|) ; with the U — 
propagator therein reducing trivially to g(u) = g(0) = — i/A in the Kondo limit. 
Using equations (5.18,21) the leading asymptotic behaviour of S(cj + ; X) (applicable for 
\uj'\ ^> 1 and all X', or X' ^> |1 + u>'\ and all uj') is thereby found to be: 



A- l T}(u + ;T) 



T) ~ — sgn(u;) tanh ( ^— ] — In 



16 

3^2 



In 2 



([vrXf +K| 2 ) 



([7rX'] 2 + |^'| 2 )^ (5.27a 



(5.276) 



For \u'\ ^> max(l,X'), the X-dependence of £(co> + ;X) is irrelevant, and 
equations (5.27) reduce to Aq 1 E r (o; + ;X) 



-sgn(u;)|fm|u/| and A E I (o; + ;X) ~ 
j^r In 2 \uj'\. Again, these coincide with the X = result, also found in || to give good 
agreement with X = NRG calculations. For X' 3> |l + u/| by contrast X' rather than u' 
controls the logarithms, which incipient divergence then dominates the imaginary part 
of the single self-energy, H 1 (uj + ; X) ~ ln 2 (X'). 
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5.2. Scaling spectra: all scales. 

The asymptotic behaviour obtained above encompasses high frequencies for all 
temperatures, as well as high temperatures for any frequency; and is independent 
of the detailed form of the transverse spin spectrum x + ~i UJ )- To obtain the LMA 
scaling spectrum on all energy/temperature scales, and in particular the low-(u;,T) 
behaviour, requires by contrast a full specification of x + ~{ UJ )- This we now turn to, 
considering as in || two related variants of the function F(u) that determines x + ~( UJ ) 
via equation (5.3). 

In the first, referred to from now on as the LMA(RPA), is obtained 

explicitly from the RPA-like p-h ladder sum equation fl2.11|) . F(u) then has the form 
equation (|5.4j) [|], and it is known |5|,|| that in the strong coupling limit U ^> 1, a — > 1 
and A — > such that Af(y) = S(y — 1). From equation (5.3a), ^x + ~( UJ ) — 8{w ~ ^m) 
thus reduces to a simple ^-function centred on u> m . The Kondo limit scaling spectrum 
arising from the LMA(RPA) then follows in closed form from equations (5.1,7,12,13), 
together with equation (5.15) for H(x). It indeed recovers precisely the scaling spectra 
obtained numerically in §4 by progressively increasing U, and illustrated in figure 4. 

The limitations of the above LMA(RPA) reside in the ^-function form of x + ~( UJ ) 
- in reality % + ~(u;) will have non-zero width, reflecting a finite lifetime for the spin- 
flip excitations. To encompass this we proceed as in ||, retaining the functional form 
equation ( |5.4j ) for F(u) (which has finite width provided a ^ 1), and employing a 
high frequency cutoff oj c to render F(u) normalizable (equation (|5.3fiQ ). The width 



parameter a is then determined by requiring that the leading low-frequency behaviour 
of the imaginary part of the T = single self-energy is recovered exactly. As explained 
in || this requires A = ?[u m / A Z} 2 (with Z the T = quasiparticle weight), which 
in turn determines a uniquely for a chosen cutoff Q c . The latter is of course arbitrary 
but, as one expects, results are not sensitive to it 0; in practice, as in ||, we choose 
uj c = 10. From now on we refer to this simply as the LMA, and recall in passing that || 
the (T = 0) HWHM Kondo scale u K = 0.691^ for both it and the LMA(RPA). 

The resultant LMA scaling spectra are shown in figure 5 for the same temperatures 
considered in figure 4 for their LMA(RPA) counterparts, T/ujk = 0.1, 0.5, 1, 5 
and 10. The differences between the two are seen to be rather minor; in fact for 
T/wk > 1 the corresponding spectra are essentially coincident for all frequencies. 
Not unexpectedly, the primary differences arise at low temperatures and frequencies, 
comparison of figures 4,5 showing that for T/u>k < 1 the Kondo resonance diminishes 
less rapidly with increasing temperature in the LMA(RPA). 

To pursue this we focus on the thermal evolution of the spectrum at the Fermi 
level, D(u = 0; T); which is related to the single self-energy by 

ttA 0j D(0; T) = [1 + AjT^O; T)}' 1 (5.28a) 

with 



2 



E^O; T) = £1(0; T) + 1 \' (5.286) 
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Figure 5. LMA scaling spectra obtained as explained in text: ttAqD(lu; T) vs w/wk 
for T/luk = 0.1 (solid line), 0.5 (long dash), 1 (short dash), 5 (point-dash) and 10 
(dotted). 



given in terms of the LMA self-energies (as follows generally from equation ( |2.13| ) 
using p-h symmetry). The resultant T/c^K-dependence of ttA D(0;T) is shown in 
figure 6, from which the LMA and LMA(RPA) are indeed seen to be essentially 
coincident for T/c^k ^ lj while differing quite significantly at lower temperatures. For 
the LMA(RPA) it is readily shown using the preceding results that the leading low- 
T dependence is £*((); T) ~ 0([T/cu K } 4 ), and hence ttA £>(0;T) ~ 1 - 0([T /uj k } 4 ). 
For the LMA by contrast, the corresponding T-dependence is of the correct form [[J 
S ! (0;T) ~ 0([T/cg>k] 2 )- The leading low T (= T/u m ) and Co (= Lj/uj m ) dependence of 
the LMA E I (u; + ; T) can in fact be obtained analytically using equations (5.7,12,13) with 
equation (5.4) for F(u), together with equation (2.13) relating E(tu + ;T) to the LMA 
T la (u + ; T). The resultant S i (cj + ; T) is thereby found to be 



T) ~ A 
1 

~ 2 



~9 71-2 ^9 

3 



A Z 



2 2 
7T 

+ T 



T 



A Z 



A Z] 2 is used, and hence 7iA D(0;T) 



where A = 

Equation (5.29) differs from the known exact result 
[T/AqZ] 2 , the exact value for which is 



1 



(5.29a) 
(5.296) 
■[T/A Z] 2 . 



only in the coefficient of 
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Figure 6. Fermi level spectrum ttAoD(uj = 0; T) vs T/luk for the LMA (solid line) 
and LMA(RPA) (dotted); the point-dash line shows the high temperature asymptotic 
behaviour equation ( 5.25| ). NRG results |2lj] for the Kondo model are also shown 
(dashed line). 



Figure 6 also shows the explicit high temperature asymptotic result for 7rAo-D(0; T), 
equation ( |5.25| ), which is essentially perfect quantitatively for T/u K > 10, and accurate 
to within 10% down to T/ujk — 5. We believe in fact that the leading high-temperature 
form from equation ( |5.25| ), txAqD(uj]T) ~ 37r 2 /[16 \t?{T/ujk)] is exact; for essentially 
the same reasons that the present theory captures the exact high temperature behaviour 
of the resistivity, as discussed in the following section. 

Finally, the LMA results for ttA D(0;T) are also compared in figure 6 to recent 
NRG calculations |21J performed directly on the Kondo model (with phost(0)J = 0.19); 
the NRG results for -D(0; T)/D(0; 0) are shown, and allow for the slight offset in the 
NRG A .D(0; 0) below its exact value of l/ir. Agreement is qualitatively good although, 
as expected from equation (5.296) where the coefficient of the leading low-T term is | of 
its exact value, the LMA ttA D(0; T) is higher than the NRG data at low-temperatures; 
and this behaviour persists over the temperature range for which NRG data has been 
published. No great leap of faith is however required to anticipate that NRG results 
for higher temperatures are likely to concur with the high-temperature behaviour of the 
present theory. 
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Figure 7. Resistivity p(T)/p(0) vs T/ujk for the Kondo scaling limit; LMA (solid 
line), LMA(RPA) (dashed). Inset: Approach to universal scaling with increasing U, 
p(T)/p(0) vs T/ujk for the LMA(RPA) with U = 4 (dotted line), 6 (dashed) and 8 
(solid line). 



6. Resistivity. 

One reason why a knowledge of the single-particle spectrum D(u;T) is important, is 
that transport properties may be obtained directly from it, as mentioned in §2. We 
consider here the resistivity p(T) (equation (|2.2j)), focussing naturally on the Kondo 
scaling regime; and for which the present theory is found to be asymptotically exact 
at high temperatures, while also recovering the correct Fermi liquid form p(T) / p(0) ~ 
1 — P[T/uk} 2 as T — > and yielding rather good agreement with NRG calculations 

flic 



First, we consider briefly the approach to the Kondo scaling limit with increasing 
interaction strength U. The inset to figure 7 shows the resultant p(T)/p(0) vs T/ujk 
obtained from the LMA(RPA), for three different U: 4, 6 and 8. From this the rapid 
approach to the limiting scaling form is self-evident. Deviations from scaling behaviour 
naturally set in for non- universal temperatures T on the order of A . This leads to 
the characteristic [/-dependent 'upturn' in p(T) evident in figure 7 (inset); which with 
increasing U moves rapidly to exponentially large vales of T/uk (wk oc exp(— 7r[//8A )), 
such that more and more of the Kondo scaling form is recovered. The main part of 
figure 7 shows the T/c^K-dependence of the resistivity in the Kondo limit, obtained 
using the analytical results of §5. Results for both the LMA(RPA) and LMA are shown 
and, as expected from the corresponding discussion of spectral asymptotics (§5.1), are 
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Figure 8. p(T)/p(0) vs T/uk for the LMA (solid line) compared to the Hamann 
approximation equation (|6.2|) (dashed line). 



coincident in practice for T/uj^ ^ 1- 

We thus consider first the high temperature behaviour of p(T). This is controlled 
by the high temperature form of D(u;T), the leading asymptotic behaviour of which 
is given explicitly by equation (|5.26|) . Inserting the latter into equation fl2.2|) for p(T), 
and transforming the integration variable therein from u to x = ou/T = oo'/T', yields 
straightforwardly the leading asymptotic behaviour for T/u^ ^> 1, 

P(T) 3vr 2 1 

p(0) ~ 16 ln 2 (T/c K ) 1 ' ) 

(where cur oc u' m is used). This is indeed the exact high-temperature asymptote 
for the S — | Kondo/s-d model, first obtained as the leading logarithmic sum of 
parquet diagrams by Abrikosov | I7|| . Note also that this behaviour mirrors closely 
the leading high-frequency asymptotics of the T = single-particle spectrum, given from 
equation flOg ) by nA D(uj;0) ~ 3tt 2 /[16 ln 2 (\cu\/tu K )]. 

Resummation of parquet diagrams for the Kondo/s-d model leads further to the 
well known Hamann approximation for p(T) [|T^] , 



P(T) 
P(0) 




ln(T/T KH ) 



[ln 2 (T/T KH ) + ^] T 



(6-2) 



which correctly recovers equation (|6.1|) to leading order for T/uk 3> 1 (Tkh oc ujg)- A 
fit of equation (|6.2j ) to our results is shown in figure 8 (with T kh /uk = 0.76). From 
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this, as found also in NRG calculations ( [IIBL see also figure 9 below), it is seen that the 
Hamann result accounts well for p(T) down to temperatures T/cjk ~ 1; while the leading 
high-temperature behaviour equation ( |6.1|) is by contrast accurate only for T/uk > 10 2 . 
We note too the parallel between the Hamann results for p(T), and equation ( |5.23| ) for 
the high-frequency behaviour of the single-particle spectrum for T = 0; where the latter 
in practice agrees well with NRG results for D(u\ 0) down to uj /uj k ~ 1 ( JT5|], see e.g. 
figure 2 therein), even though its leading high-frequency asymptote oc l/ln 2 (|u;|/u;K) is 
accurate only for u/uk > 10 2 . 

We turn now to the low-temperature asymptotics of p(T), which as first shown by 
Nozieres [ 19| has the characteristic Fermi liquid form 



P{T) 1 
P(0) 



T 



A n Z 



(6.3) 



with the exact coefficient b = tt 2 ; and where AqZ (with Z the T = quasiparticle 
weight) is related to the conventional Kondo temperature T K , defined such that the static 
impurity susceptibility Xim P = (<?Pb) 2 /4T k , by T K = jA Z. The T 2 -expansion of p(T) 
may be extended to higher order via the methods of boundary conformal field theory ||23|| , 
and in principle by the more generally applicable renormalized perturbation theory 
|24],[25|]. The former has been carried out exactly up to (T/AqZ) 6 [^3|; but this impressive 
feat barely extends the applicability of the T 2 -expansion beyond T/A Z ~ 1CT 1 , up to 
which the Fermi liquid result equation ( |6.3|) is accurate in practice. 

As shown in the leading low-T behaviour of p(T) for the AIM generally is 
related to the low-(u;, T) behaviour of the single-particle spectrum by 

P(T) . J\ , W fd 2 A(u;0)\ | T 2 ( d 2 A(0;T) \ j (g4) 



P(0) I 6 V ) W=Q 2 V dT* /T=QJ 

where A{uj\ T) = 7tAqD(uj;T). In the Kondo limit the exact results are P, |16[ 
A(cu;0) ~ 1 - 1{lu/A Z] 2 and .4.(0; T) ~ 1 - ^-[T/A Z] 2 ; equation Q thus follows. 
The leading low-T behaviour of p(T) from the present local moment approach can 
likewise be determined. It recovers the Fermi liquid form equation ( |6.3j ), scaling in 
terms of T/A Z, although the coefficient b is not exact. For the LMA, A(u;0) to 
leading order (and hence the first non-trivial term in equation ( |6.4| )) is exact, while 
.4.(0; T) = [1 + Aq^^T)]- 1 ~ 1 - ^-[T/A Z] 2 (from equation QQ91D ); so the 
coefficient b = 2vr 2 /3. And for the LMA(RPA), where A(u; 0) ~ 1 - [u/A Z] 2 and 
A(0;T) ~ 1 - O([T/A Z} 4 ) (§5.2), b = vr 2 /3. The difference between the LMA(RPA) 
and LMA at low temperatures is of course seen directly in figure 7. 

In figure 9 the LMA p(T)/p(0) vs T/T K is compared directly to the NRG results of 
Costi, Hewson and Zlatic |nj. The latter were obtained for the AIM with [7 = 4; and 
the final two/three high-T NRG points shown correspond to T ~ O(Aq), thus exhibiting 
the characteristic non-universal upturn from Kondo scaling behaviour mentioned above 
in relation to figure 7 (inset). We regard the agreement with NRG as rather good, 
the more so given the relative simplicity of the local moment approach and its ready 
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Figure 9. p(T)/p(Q) vs T/T K (with T K defined such that Xim P = (g Mb) 2 / '4T K ) . LMA 
(solid line) and NRG results (l|] for the AIM with U = 4 (diamonds). The exact 
high-T asymptote p(T)/p(0) ~ 3tt 2 /[16 ln(T/T K )]fij§ is also shown (dashed line); 
as are results from a recent low-T approximation |2q| exploiting integrability (dotted 
line). See text for discussion. 

extendibility to a much wider range of problems than the integrable AIM. Figure 9 also 
shows results from a recent approach [fj6| that, while unable to access single-particle 
dynamics, focusses directly on transport (strictly the linear differential conductance, 
which is however known to differ negligibly from p(T) pip. This method exploits the 
integrability of the AIM via an approximate treatment of the Bethe ansatz equations, 
and is designed to capture the low-T behaviour. That is does so very well up to T/T K ~ 1 
is evident from figure 9; although it does not appear to recover the high-temperature 
Abrikosov asymptote equation ( |6.1| ) (referred to in |26| as the 1-loop RG result after its 
recent essential rediscovery in |27j via an elegant RG approach). The latter (specifically 
37r 2 /[16 ln 2 (T/T K )]) is also shown in figure 9 and is recovered by the present work, in 
practice for T/Tk > 10 2 . 

7. Conclusion. 

The subject of this paper has been single-particle dynamics of the Anderson impurity 
model, pursued via the local moment approach extended to finite temperature, and with 
a natural emphasis on the strongly correlated Kondo scaling regime. The approach 
yields thereby a rich, and rather successful description of the thermal destruction of 
the Kondo resonance and hence of transport properties such as the resistivity. This 
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augurs well for the future, since the LMA is neither technically difficult nor confined 
to the integrable, metallic AIM considered here; suggesting its practical viability as a 
potentially powerful route to dynamics and transport properties of a wider range of 
quantum impurity models, as well as e.g. the periodic Anderson and Hubbard [pEf 
models within the framework of dynamical mean field theory. Problems of this ilk are 
currently under study, and will be considered in subsequent publications. 
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